Dynamics of correlations in atomic Bose-Einstein condensates 
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The Gross-Pitaevskii equation has been extremely successful in the theory of weakly-interacting Bose- 
Einstein condensates. However, present-day experiments reach beyond the regime of its validity due to the 
significant role of correlations. We review a method of tackling the dynamics of correlations in Bose condensed 
gases, in terms of noncommutative cumulants. This new approach has a wide range of applicability in the areas 
of current interest, e.g. the production of molecules and the manipulation of interactions in condensates. It also 
offers an interesting perspective on the classical-field methods for partly condensed Bose gases. 

Eight years after the first experimental achievement of Bose-Einstein condensation of atoms fjj, the field of quantum- 
degenerate gases is not losing its pace. On the contrary, it continues to open new andpromising avenues for research, with 
the latest hot topics including quantum vortices OyL Bose condensation of molecules |4] and prospects of elucidating pairing 
mechanisms with the help of ultracold Fermi gases |5]. 

Given the enormous potential and diversity of the field, it comes as a surprise that most of the experiments performed so 
far with Bose-Einstein condensates have been modelled almost perfectly by mean-field theory. The famous Gross-Pitaevskii 
^sO ' equation |6, 7] has proved to be a genuine workhorse in the area, being employed in a myriad of contexts |8]. An interesting and 
relatively new addition to the plethora of applications is a classical-field method for the dynamics of interacting Bose gases at 

non-zero temperatures |3|l3|ni[l3li3IIllll. 

Only recently have new types of experiments been performed which necessitate the use of more advanced theoretical tech- 
ly-j ■ niques. A striking example is provided by an experimental demonstration of a quantum phase transition from a superfiuid to a 
y^Q ' Mott insulator with bosons in an optical lattice II 611. In this case, the mean field vanishes completely in the insulating phase. 
Mean-field theory also fails in a strongly-interacting regime accessible through magnetically tunable two-body scattering in 
the vicinity of a Feshbach resonance 1171 Hal . Very recentl y, a number of experiments have explored the association of ultra- 
cold molecules through crossing a Feshbach resonance 1 191 1201 12111 . Again, in this situation mean field approaches provide an 
1 approximate description of the processes involved in limiting cases only j2*3l . 

Apart from the availability of new experimental results challenging the applicability of the lowest order, mean-field techniques, 

S a very interesting general issue is to what extent the time-dependent Gross-Pitaevskii equation describes a dynamical response 
of a condensed gas to perturbations. Experiments operating in the vicinity of a Feshbach resonance provide just such an example 
where external perturbations are violent and may lead to an almost complete depletion of the condensate. However, one can 
imagine other, more general situations, e.g. a rapid variation of a trapping potential, where losses can also be expected to be 
O ■ significant. 

In this paper, we review a recently developed method capable of handling the dynamics of a Bose condensed gas beyond mean 
J> , field theory in a self-consistent, number-conserving way 1 23] . The method is based on the expansion of correlation functions 
in terms of noncommutative cumulants 12411 . After sketching the general ideas behind this new approach, we describe a test 
of its predictions against the experimental data in a particular case of a Feshbach resonance crossing performed with a 85 Rb 
condensate at JILA l25il . Furthermore, we provide an outlook for other possible applications of the cumulant approach. 

We start from a very general form of a many-body Hamiltonian in its second-quantised form: 
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H = J dx ^( x)i j 1B (x)^( x ) + - / dxdy ^ t (x)V^ t (y)^(x - y)V(x)V(y) , (1) 

where ?/>(x) is a field operator for atoms, obeying the bosonic commutation relations: [i/j(x.), V> (y)] — <5( x — y) and 
[ip(x),ip(y)] = 0. V(x — y) is the two-body interaction potential. The one-body Hamiltonian, which may contain an ex- 
ternal (e.g. trapping) potential, is defined as follows: 



H 1B = -h 2 V 2 /2m + V^t ■ (2) 

All physical properties of a many-body system, e.g. a gas, are determined by correlation functions, i.e. expectation values of 
normal ordered products of field operators in the state at time t (denoted hereafter as (...)*) |26]. A usual way of proceeding 
with the problem of the time evolution of correlations consists in writing the hierarchy of the equations of motion, coupled due 
to the quartic interaction term in the Hamiltonian. Here we describe a different procedure, which involves noncommutative 
cumulants. Given a set of bosonic field operators B\, B 2 , B 3 , ... the cumulants (denoted as (. . .) c ) are defined recursively in the 
following way: 
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(Bi)t = (Bi) c , 
{B 1 B 2 ) t = (B 1 B 2 ) C + (B 1 ) C (B 2 ) C , 
{B 1 B 2 B 3 ) t = (B 1 B 2 B 3 ) C + (B 1 ) C {B 2 B 3 ) C + (B 2 ) C (B 1 B 3 ) C + (B 3 ) C (B 1 B 2 ) C + (B 1 ) C (B 2 ) C {B 3 ) C , 

: (3) 

It follows from the above definition that the cumulants capture the essential correlations in the system by subtracting factoris- 
able contributions from the correlation functions. They decrease with an increasing order, provided that the system is not too 
far away from the interaction-free equilibrium, which is typically not true for correlation functions. In particular, for an ideal 
gas in thermal equilibrium all cumulants containing more than two field operators vanish, in agreement with Wick's theorem of 
statistical mechanics. Thus, higher-order cumulants measure the deviation of the system from the interaction-free situation. 

Specifying the cumulant expansion in the case of a partly condensed gas of bosons, one expects a special role to be played by 
the following three lowest-order cumulants: 



*(x,i) = (^(x))t, (4) 
*(x,y,t) - ^(y)^(x)) t -*(y,i)*(x,*), (5) 

r(x, y ,t) = (^ t (y)^(x)) i -**(y,t)*(x,i). (6) 

The first-order cumulant ^(x, t) is the mean field, while the second order cumulants $(x, y, t) and T(x, y, t) will be referred 
to as the pair function and the density matrix of the non-condensed fraction, respectively. We note that the commonly used 
Gross-Pitaevskii formalism with the contact potential would involve factorising all correlation functions, which corresponds to 
neglecting all cumulants of order higher than one. 

The next step is to write the equations of motion for the cumulants. In order to handle the resulting infinite hierarchy of 
them, we define the following truncation scheme of order n. The equations of motion for cumulants containing not more than 
n operators are kept exactly. In the equations of motion for the cumulants of order n + 1 and n + 2 we include free evolution 
only, i.e. we neglect those products of normal-ordered cumulants that contain n + 3 and n + 4 field operators. As a result, the 
equations of motion for the cumulants of order n + 1 and n + 2 can be solved formally and the solutions can be substituted 
back to the equations of motion for cumulants of order not greater than n. The consequence of such a procedure is that bare 
microscopic potentials are connected to transition matrices in all equations of motion up to order n. 

In what follows, we present the first-order (n — 1) cumulant approach for bosons |23|. In this case the relevant equations of 
motion, after the first-order truncation, take the following form: 

aJU( x ,<) = Jf 1B (x)*(x,t)+ j dyy(x-y)[(^( y )^( y )^(x)) c + *(y,t)r(x,y,i)+vl/(x,<)r(y,y,t)] 

+ y"dyy(x-y)**(y,t)[*(x,y,t)+*(x,t)*(y,t)] , (7) 

iftJU(x,y,i) = fl- 2 B(x,y)*(x,y,t)+K(x-y)*(x,t)*(y,*), (8) 

where H 2 b{x, y) = -Hib(x) + i?iB(y) + V(x~ y). The first-order approach also includes the appropriately truncated equations 
of motion for T(x, y, t) and {ip^ (x)ip(y)'tp(z)) c . We note |23], however, that if T(x, y, t) and (ip^ (x.)ip(y)ip(z)) c vanish initially, 
they will not evolve. Thus, we neglect these two cumulants on the right-hand side of Eq.Q. Back substitution of the formal 
solution of the linear equation of motion for $(x, y, t) leads to a closed equation of motion for the mean field ^(x, t): 

O ft 

i%— *(x,t) = Hi B (x)*(x,t)+**(x,t) ^ dti(2 7 r?i) 3 (0|T 2 ( + ) (t,i 1 )|0)* 2 (x,t 1 ), (9) 



where 



T^\t, h) = V6(t - h) + VG { 2 1 ] (t, t x )V (10) 



is a retarded two-body time-dependent transition matrix. The time-dependent transition matrix is expressed in terms of the 
retarded two-body Green's function G 2 b' (*> *i)> which is defined in a usual way through 



ih^- t -H 2B )G ( 2 l ) (t,t 1 )=5(t-t 1 ). (11) 
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The two-body physics enters the first-order cumulant approach through the two-body transition matrix. Thus, the binary colli- 
sions are described non-perturbatively. In Eq.(|5), the transition matrix is sandwiched between two plane waves of zero (relative) 
momentum (denoted as |0)). This reflects the short range of the microscopic interatomic potential V as compared to the length 
scales relevant to the mean field ^(x, t). In practical terms, Eg. (II 1> . describing the two-body dynamics, can be solved by making 
use of a separable form of the two-body interaction potential (for details, see II221I27II ). which is a rather good approximation for 
low-energy collisions. In the end, the two-body problem is described in terms of a few collisional parameters (e.g. the scattering 
length or the scattering resonance width), which can be determined experimentally or theoretically. 

Eq.(|9} has the form of a non-Markovian non-linear Schrodinger equation. It can be shown that in the Markov approximation 
it reduces to the well known time-dependent Gross-Pi taevskii equation J23ll . We note in passing the limitations of the above for- 
mulation of the first-order cumulant approach. Neglecting the initial one-body density matrix of the non-condensate T(x, y, to) 
implies that the initial state must be described by a mean field. The above derivation also assumes that the initial pair function 
$(x, y, to) vanishes, so no account of collisions before the initial time is taken. Both assumptions are well satisfied in most 
present-day experiments. 

Although r(x, y, t) does not appear explicitly in the equation of motion (|9j (which reflects the assumption that the non- 
condensed fraction does not affect the dynamics of the mean field), the atoms lost from the condensate are transferred into the 
non-condensate. Having solved for the evolution of the mean field and the pair function, we can now use them to determine 
the evolution of the non-condensed fraction by retaining only these terms in the full equation of motion for T(x, y, t) ¥Fh . The 
resulting dynamical equation has the following form: 

^r(x,y,t) = |if 1B (x)r(x,y,f) + J rfzy(x-z)$*(z,y,t)[$(x,z,i)+*(x,<)*(z,t)]|-|x^y| . (12) 

It can be shown that for an initially undepleted condensate, i.e. r(x, y, to) — 0, the solution of Eq. J12t can be expressed in 
terms of the pair function as 



r(x, y, t) = J dz $(x, z, <)$* (y, z, t) . (13) 

Consequently, the non-condensed fraction stems from the build-up of pair correlations. In turn, these correlations can appear 
either in the scattering continuum (as correlated pairs of relatively fast atoms) or in the bound part of the two-body spectrum, as 
molecules |27]. Given Ea.(ll3>. the total number of atoms N can be shown to be a constant of motion: 



N = J dx[r(x,x,i) + |*(x,t)| 2 ] . (14) 

As an application of the cumulant method we study the experiment of Ref. 12511 performed with a 85 Rb condensate at JILA. 
In this experiment, the magnetic field was swept across a Feshbach resonance, starting from the side of the positive scattering 
length. A loss of atoms from the condensate was observed, which depended on the ramp speed of the magnetic field strength. In 
applying the cumulant approach to this case, special care has to be taken to describe properly the underlying two-body dynamics 
in the vicinity of a Feshbach resonance 1 28] . We have solved the non-Markovian non-linear Schrodinger equation © to simulate 
the actual experimental procedure. At the end of the ramp we record the decrease in the norm of the mean field \&(x, t) and 
compare it with the experimentally measured condensate loss. The results are shown in Figure^ 

We have employed the non-Markovian equation of the first-order cumulant approach to study other situations where significant 
loss of condensate atoms due to a dynamic perturbation is observed. The method has been successful in predicting the rate of 
atom loss in four-wave mixing experiments |23|. We have also analysed 1 27] the experiment pioneering the use of magnetic-field 
pulses for the association of molecules, where coherent atom-molecule oscillations have been observed 1 18]. In this context, the 
physical consequences of an extremely long-range nature of molecules created have been pointed out 1 29] and the process of 
adiabatic association of dimer molecules through Feshbach resonance crossings has been analysed in detail 1 22] . The second- 
order approach has been used to determine the mean field energy associated with three-body collisions in Bose condensates 

There are other situations of experimental and theoretical interest where the mean-field dynamical description can be expected 
to fail due to the significance of pair correlations. In the context of optical lattices, the cumulant approach could be a useful 
tool in describing the dynamics of the superfluid bosonic gas in the vicinity of the quantum phase transition, where quantum 
depletion starts to grow. An area relatively unexplored experimentally is the dynamic manipulation of magnetically-tunable 
interactions, not necessarily in the close vicinity of a Feshbach resonance 1 3 1 ] . The question of adiabaticity of the evolution with 
respect to particle loss from the condensate is an important one in this context. Finally, one can study situations where a Bose 
condensate is subjected to rapid variations of an external (e.g. trapping) potential. One of the interesting proposals in this field 
is a delta-kicked oscillator setup, employed in the studies of quantum chaos H32I1 . 



4 



80 



5? 



■60 



o 
-J 



o 
U 



• expenment 
■ theory 



50 100 150 

Inverse Ramp Speed [us/Gauss] 



200 



FIG. 1: Loss of atoms from a 85 Rb condensate at the end of a linear magnetic field ramp across the 155 Gauss Feshbach resonance as a 
function of the inverse ramp speed. The circles with error bars are the experimental data from Ref. |25]. The squares are determined from 
Eq. |5J using a cylindrical trap potential with the experimental trap frequencies. 



From the viewpoint of the classical field methods, the cumulant approach also may offer an interesting perspective. The non- 
Markovian equation can be viewed as a consistent way of including quantum corrections in the classical-field treatment. As 
such, it can be used to test the validi ty of rep lacing quantum field operators by complex-valued functions, which is the backbone 
of the latter approach BE [l(J E III El • 

To conclude: in this paper we have reviewed a new method of tackling the dynamics of correlations in a Bose condensate, in 
terms of noncommutative cumulants. This framework delivers an accurate but computationally feasible approach to the beyond- 
mean-field evolution in a Bose-condensed gas. The method has been used to address several issues of current interest. Finally, 
we have pointed out several new applications and a possible link to the classical-field methods for partly condensed Bose gases. 
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Research Fellowship of the Royal Society (T.K.), Deutsche Forschungsgemeinschaft (T.G.). K.B. thanks the Royal Society and 
the Wolfson Foundation. 
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